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We study the effect of electrostatic interactions on the distribution function of the end-to-end 
distance of a single polyelectrolyte chain in the rod-like limit. The extent to which the radial 
distribution function of a polyelectrolyte is reproduced by that of a wormlike chain with an adjusted 
effective persistence length is investigated. Strong evidence is found for a universal scaling formula 
connecting the effective persistence length of a polyelectrolyte with the strength of the electrostatic 
interaction and the Debye screening length. 



PACS numbers: 82.35.Rs, 87.15.La, 36.20.-r, 82.35.Lr 



I. CONFORMATIONAL STATISTICS OF A 
POLYELECTROLYTE IN THE ROD-LIKE LIMIT: 
BACKGROUND AND GENERAL 
CONSIDERATIONS 

Because of their fundamental importance in many key 
biological processes, electrostatically charged polymers, 
or polyelectrolytes (PE's), have been the subject of in- 
tense research recently. These polymers often have stiff 
structures due to their Coulomb self-repulsion, and this 
structural property lies at the heart of their biological 
functionality. For example, the cytoskeletal network of 
actin filaments plays a crucial role in the recovery of eu- 
karyotic cell shape in the face of the stresses imposed by 
cell movement, growth, and division ^ ^. Further- 
more, the storage properties and accessibility of DNA 
is known to be both constrained and controlled by its 
stiffness |Q, |[ . Clearly, a complete study of the elastic 
properties of these charged biopolymers, which is now 
facilitated by single filament imaging and manipulation 
techniques, is crucial to an understanding of their role in 
nature ||^. 

In the case of neutral polymers, the wormlike chain 
(WLC) model of Kratky and Porod provides a power- 
ful and convenient characterization of flexibility through 
the persistence length, which quantifies the correlations 
of unit tangent vectors parallel to the chain. Odijk, Skol- 
nick, and Fixman (OSF), have applied this notion to PE's 
through the introduction of an "electrostatic" persistence 
length if, 0]. They derived the relationship 



tOSF 



(1) 



where is the Debye screening length (a measure of 
the ionic strength of the solvent), and /3 = i's/^^ is the 
strength of the electrostatic interaction. In the latter 



quantity £b — /^ksT is the Bjerrum length with e the 
dielectric constant of the ion-free solvent, and b the av- 
erage separation between neighboring charges along the 
PE. 

The notion of electrostatic persistence length, how- 
ever, is believed to have shortcomings, as it replaces the 
many length scales in a PE by a single effective persis- 
tence length. The first note of caution is due to Odijk 
himself, who pointed out that if one tries to make use 
of the concept of an effective WLC approximation for 
a PE of finite length L, the corresponding persistence 
length should be X-dependent [||. Barrat and Joanny 
later noted that, in fact, each deformation mode of a PE 
with a given wavelength has a different effective persis- 
tence length Moreover, the wavevector dependence 
of the rigidity causes deviations from the linear depen- 
dency 4 oc /3 |10[ |l). 

These observations — which point to the complications 
that may arise if one insists on an effective WLC scheme 
for the characterization of PE's — follow from the study of 
the averages {R) and (i?^), where R ~ r(L) — r(0) is the 
end-to-end distance of the PE. It is not at all clear that a 
description based on the matching of first or second mo- 
ments of the end-to-end length distribution is equivalent 
to a comprehensive study of conformational properties 
of a PE. An appealing alternative is a study of the full 
radial distribution function of the end-to-end distance of 
PE's, which provides an excellent gauge of the conforma- 
tion statistics of polymers in general [|l2j and stiff chains 
in particular ijl^ . This quantity, which can be measured 
experimentally through fiuorescence microscopy, allows 
one to determine whether or not the WLC model accu- 
rately reflects the properties of the polymer segments. If 
the model has been shown to apply, then the fit of the 
distribution to the data directly yields the persistence 
length and, hence, the elastic modulus of the polymer. 
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In this article, we calculate the end-to-end distribution 
function of a charged inextensible chain in its rod-like 
limit. This distribution is then utilized to investigate the 
validity of a WLC with an adjusted persistence length 
as a model for a PE. We find that when the Coulomb 
interaction is only a perturbation to the mechanical stiff- 
ness of the chain, or when Debye screening is sufficiently 
strong, the full radial distribution of the PE can be col- 
lapsed onto that of a WLC. However, in the case of sizable 
Coulomb interactions and when the screening length is 
insufficiently short, we find that the distribution cannot 
be collapsed onto that of any WLC. To quantify the de- 
viations from WLC behavior, we define an effective £p 
as the persistence length of the WLC that is most likely 
to have the same end-to-end distance as the PE chain. 
We discover strong evidence of a universal scaling func- 
tion relating the electrostatic persistence length, i^, to 
the intrinsic persistence length, £pQ, the Debye screening 
length, K and the strength of electrostatic interaction, f3. 

Because of the inextensibility of the PE's under con- 
sideration we adopt the Kratky and Porod WLC 
model to describe the bending energy of the chain. In 
this model, a polymer is represented by a space curve 
r(s) as a function of the arc length parameter s. The to- 
tal energy of the chain, which is the sum of the intrinsic 
elasticity and the electrostatic energy can be written as 



H 

knT 



ds 



dsds' ■ 



-K|r(s)-r(s')| 



r(s)-r(sOI 



(2) 



where t(s) is the unit tangent vector, and ipo is the in- 
trinsic persistence length of the PE. We do not take into 
account the fluctuations in the charges localized to the 
chain and in the counterion system that can give rise to 
attractive interactions leading to chain collapse 1 14 . Be- 



cause the chain is in its rod-like limit, excluded volume 
does not play a role. The double integral on the right 
hand side of Eq. is cut off when |s — s'| < b, where 
b is the spacing between charges on the PE. The quan- 
tity b also represents the intrinsic "coarse-graining" of 
the effective Hamiltonian, in that structure on a smaller 
length scale is not encompassed by the model embodied 
in Eq. (^. In addition, the value of b represents a lower 
limit on the magnitude of the intrinsic persistence length 
£pO, in that there is no physical meaning to a persistence 
length that is exceeded by the smallest length scale in 
the system. 

The end-to-end distribution function is defined as fol- 
lows: 

G(r) = (5(r-R)), (3) 

where R — r(L) — i"(0). The average in (||) is over an 
ensemble of PE chains. The function G(r) is, then, the 
probability that a given chain in the ensemble will have 
an end-to-end distance equal to r. We make use of the 
procedure Wilhelm and Frey have implemented to calcu- 
late the end-to-end distribution function for inextensible 



neutral polymers We focus on Inextensible chains 

with ipo L or £pQ ^ L. In both cases, our study 
is limited to regions where the combination of intrin- 
sic stiffness and repulsive strength of the Coulomb in- 
teraction keeps the chain in its rod- like limit That 
is, the combination of intrinsic persistence length and 
Coulomb repulsion acts to keep the chain nearly straight 
over its entire length, so that its end-to-end distance 
does not greatly deviate from its total arc length. We 
parameterize the contour in terms of the tangent field: 

t(s) = {ax{s), ay{s), 1)1 ^1 + a1{s) + a2(s), and then re- 
tain only terms up to second order in the a's in the mea- 
sure factor, TJ , and in the argument of the delta-function 
m Eq. (|), which is rewritten in terms of its Fourier rep- 
resentation. Expanding ax{s) and ay{s) in a cosine series, 
as mandated by the open-end boundary conditions, and 
making use of the relationship 

r(s) -r(s') = / dxt{x), (4) 

we obtain 



in which is a normalization constant. The A„'s are 
eigenvalues of the matrix 



T= (n^)2(£poA)I + /3LE, 



(6) 



where I is the unity matrix and E is the electrostatic 
energy matrix in the cosine basis set. Contour integration 
then yields 



G(r)=iV^/(n)e- 



A„(l-r/L) 



where 



/(n) = n 1/(A. 



Xn). 



(7) 



(8) 



The large eigenvalues of the matrix T are dominated 
by the diagonal terms {mr)"^ {£po / L) . In other words, the 
effect of the electrostatic interaction is swamped by semi- 
flexible energetics at short length scales. In the case of 
a neutral polymer with an intrinsic persistence length 
that is comparable with the total length of the polymer, 
the above series converges very rapidly, as noted in Ref. 
Jl3| . When the chain is charged, and in particular, when 
the stiffness is due predominantly to electrostatic effects, 
more terms in the series must be preserved in order to 
obtain a stable answer for G{r). In the calculations de- 
scribed here, we truncate the matrix T at a size much 
greater than the number of terms needed to obtain an 
accurate answer for the series of Eq. (^). This is be- 
cause a larger dimension of the truncated matrix leads 
to more accurate values for the lower eigenvalues (which 
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participate in the sum), and a more accurate f{n). In 
practice, we increased the dimension of the matrices un- 
til the effect (on the end-to-end distribution function) of 
a further increase in the size of T by a factor of two was 
less than a part in a thousand. A restriction on the size 
of the matrix arises from the requirement that the coarse 
graining length b not exceed the smallest wavelengths ap- 
pearing in the cosine basis set. If the length of the PE 
is L, this means that the size, iV, of the basis set satis- 
fies N < L/b. At no point in our calculations was this 
inequality violated. 

We set b/L = 10^"^. In this case, the distribution func- 
tion of PE's depends on three independent dimensionless 
parameters £po/L, (3L, and kL. Figure |l| illustrates the 
effect of the screened electrostatic interaction on the dis- 
tribution function of a chain with ipo/L = 0.5 at different 
salt concentrations and, hence, different values of k. As 
illustrated in the figure, the end-to-end distance of the 
chain on average becomes shorter upon a decrease of the 
screening length (for fixed PL). The quantity PL, which 
depends on the solvent dielectric constant and the charge 
density on the PE chain, is equal to 100 for all the three 
curves in the figure. 

We now seek to determine under what conditions the 
end-to-end distribution function of a rod-like PE is satis- 
factorily reproduced by that of a WLC with an adjusted 
persistence length. We also scrutinize the adjusted per- 
sistence length of a WLC whose end-to-end distribution 
is the closest fit to that of a rod-like PE to determine 
whether or not this persistence length is in accord with 
existing predictions for the dependence of the electro- 
static persistence length on properties of the PE chain. 
In the cases we consider here the combination of intrinsic 
stiffness — as parameterized by Ipo/L — and of charging — 
quantified by PL — provides enough stiffening of the PE 
to guarantee rod-like behavior in all length scales. 



II. REGIMES IN WHICH THE PE BEHAVES 
LIKE A WLC 



These regions are associated by weak charging of chain 
for small values of kL ( kL ~ 1 ). As we increase 
screening, kL, the collapse of the two distributions can 
be achieved for higher charge density. Figure ^ illus- 
trates two examples of this regime. In plot (a), the De- 
bye screening length is short, while PL, the strength of 
coupling, is high. In plot (b), the Debye screening length 
is relatively long, but the electrostatic coupling is very 
weak. In both cases, we find that the distribution func- 
tion of the PE's collapse onto WLC's with adjusted per- 
sistence lengths that follows Odijk and OSF predictions 
(Eqs. §) and d)) . 

Our calculations verify that under physiological condi- 
tions {k = 1 nm~^), the distribution functions for rod- 
like DNA segments {L^lOO nm) as well as those of stiffer 
actin filaments also collapse onto the end-to-end distribu- 
tion for neutral chains with an effective persistence length 
given by Eq. (|l|). 

Whenever there is a virtually perfect collapse of the 
distribution function of a PE onto that of a neutral chain, 
the persistence length of the neutral chain follows Odijk's 
prediction, in that, £p = +^pO: where £p is the effective 
persistence length of the charged chain, and £e ~ ^odijk 



tOdijk 
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(kL)2 {kLY 



{kL)2 [nLf 



(9) 



which reduces to iosF = P/^k? for large kL It is 

important to note that the expression for -^odijk |p| was 
derived under the assumption that the contour length of 
the chain, L, is of the order of its intrinsic persistence 
length, £pQ and that the electrostatic interaction has the 
limited effect of only "perturbing" the WLC shape of the 
chain Our results confirm the validity of OSF and 
Odijk formulas in the regime in which they are expected 
to be correct. 



We find that there exist regimes in which the confor- 
mational statistics of a PE chain in the rod-like limit are 
identical to those of a WLC with an adjusted persistence 
length. For such cases, PE and WLC distributions are 
indistinguishable to the naked eye. In these regimes, the 
persistence length is as predicted by Eq. (|l|), or, when kL 
is not significantly larger than unity, by Eq. (^ . We are 
able to divide this regime into two different categories. 

A. Intrinsically stiff chains {Ipo ~ L) 

Whenever the electrostatic interaction plays a pertur- 
bative role in the chain energetics, it is possible to obtain 
a satisfactory collapse of the PE distribution onto that of 
a WLC. One noteworthy feature of this regime is that the 
electrostatic persistence length always satisfies £e < £po- 



B. Intrinsically flexible chains {£po <^ L) 

We are also able to identify regimes in which the ra- 
dial distribution function of a PE that is long compared 
to its intrinsic persistence length {£pQ/L ^1) collapses 
almost perfectly onto that of an uncharged WLC. Sub- 
stantial charging is required to enforce the rod-like limit 
for such chains. It is important to note that despite the 
substantial charging of the chain, the ratio of the length 
of the PE to the screening length, kL, must be sufficiently 
large. This requirement is essential in order to minimize 
the strong influence of the end effects on conformational 
statistics of charged chains. 

The important characteristic of this regime is that 
£pO <C £e and thus electrostatic energy no longer plays 
a perturbative role. Nevertheless, if kL is sufficiently 
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large we observe collapse of the two distributions. Fur- 
thermore, we also find that Odijk's formula works to 
near perfection in those cases. For example, in the ex- 
treme instance of a very short intrinsic persistence length, 
(■po/L — 0.01, a high degree of charging (3L = 36000, and 
a screening length that is short compared to the chain 
length, kL = 100, we obtain an end-to-end distribution 
that is nearly identical with that of a WLC with an ap- 
propriate £e- In this case, le/L = 0.894 which is in near- 
perfect agreement with ^odijk- The accuracy of Odijk's 
formula when ^ is not at all obvious, as OSF 
was derived in the regime where electrostatic interaction 
plays a perturbative role. 

We emphasize that the reason for the high quality of 
the match with OSF in this regime is different from the 
reason for the corresponding result obtained by Khokhlov 
and Khachaturian for weakly charged flexible chains. 
In our case, the chain is stiff in all length scales; the 
possibility of renormalizing the length and/or charge is 
thus excluded in our formulation. 



Our general observation is that the ratio of £e extracted 
from the effective WLC distribution to Odijk's persis- 
tence length correlates with the quality of the fit of an 
effective WLC end-to-end distribution to that of a PE. 
When this ratio is equal to one, the PE is described to 
a high degree of accuracy in terms of a WLC. As this 
ratio decreases, the deviation becomes more pronounced. 
Figure ^ displays a diagram which delineates the quality 
of the fit. The line that is used in the Figure to sepa- 
rate the two regimes corresponds to ^e/^odijk = 0.58. As 
indicated in this Figure, for fixed kL, when f3L is be- 
low a certain value the PE behaves like a WLC, while 
for larger (3L there is a substantial difference between 
the electrostatic persistence length of a rod-like PE and 
Odijk's prediction. As illustrated in Fig. |[ the deviation 
of ^e/L from ^odijk/-^, with increasing [3L is more pro- 
nounced at lower values of K.L. The Figure also shows 
that the electrostatic persistence length of a rod-like PE, 
le/L, depends on ipo/L. This dependence is not present 
in Odijk's formula. 



III. REGIMES IN WHICH THE 
CONFORMATIONAL STATISTICS OF A PE 
DIFFER FROM THOSE OF A WLC 

No matter how short the range of electrostatic inter- 
actions is as the result of Debye screening, at sufficiently 
strong charging, the PE end-to-end distribution differs 
significantly from that of a WLC. The emergence of a 
difference between the two distribution is accompanied 
by a divergence between the effective persistence length 
characterizing the conformational statistics and the pre- 
dictions of either Eq. (|l|) or (^. Figure ^, displays the 
PE end-to-end distribution (solid curves) along with the 
modified WLC distribution (dashed curves) in a case in 
which it is possible to obtain a nearly perfect fit (plot b) 
and in a case in which the best fit not nearly as good 
(plot a). In both cases fit was obtained by matching the 
location of the maxima of the two distributions, and the 
electrostatic persistence length attributed to the PE dis- 
tribution is that of the WLC associated with the dashed 
curve. The ratio of electrostatic persistence length ex- 
tracted from the distribution function to Odijk's persis- 
tence length is, -^e/^Odijk = 0.99 in the case of plot b. For 
the inset graph, plot a, in which the two distributions 
differ noticeably, this ratio is ^e/^Odijk = 0.89. 

It is important to note that a rescaling of the back- 
bone length of the PE is not an acceptable stratagem for 
improving the agreement between the PE radial distri- 
bution and that of a WLC. This is because the backbone 
length is essentially fixed by the rod-like chain condition. 
The shortening and thickening that is associated with 
intermediate blob-like structures ||T^ will not occur. In 
fact, we have observed that a reduction in the effective 
value of L actually degrades the quality of the correspon- 
dence between the conformational statistics of a PE and 
those of the corresponding WLC. 



IV. UNIVERSAL BEHAVIOR 

Our general observation of the curves of te/L versus 
f3L is that they asymptote to a power law of the form 
te oc {(3LY~^ where the exponent x{kL) lies between 
zero and one. In search of a possible universal behavior, 
which might relate ip^/L, j3L and kL to te/L, we have 
rescaled our graphs for different values of the parameters. 
The details of rescaling will be presented in a forthcom- 
ing publication The dependence on charging, (3L, 
of the ratio of the persistence length of a WLC and the 
prediction for that quantity embodied in Eq. (b|) can be 
systematized in terms of a universal formula psj. This 
formula points to the existence of a scaling mechanism. 
At this point, we are not able to provide an explana- 
tion for this mechanism or indicate a possible underlying 
basis for it. Figure ^ contains plots of the ratio of the 
electrostatic persistence length of PE distributions to the 
predictions of Odijk in Eq. (|^). As shown in the figure, 
all rescaled plots collapse on to a single curve correspond- 
ing to the following crossover formula for €e 



KOdijk 



1 + {PIPoY 



(10) 



This formula provides a remarkable fit to our data when 
the exponent x ~ 0.4, as exemplified in Fig. ^. The 
quantity (3q is an increasing function of tpQ and k. Our 
most striking result is that for all investigated values 
of kL and/or Ip^jh^ the exponent x remains fixed at 
~ 0.4, independent of all the other parameters. This 
function encompasses all the regimes described above. 
We are able to obtain the value of /3o from our data. The 
crossover discussed in the previous section occurs when 
/3 ~ 0.446/3o which corresponds to — 0.58£odijk in Eq. 
(|To|). The curve in Fig. ^ indicates the values of f3oL at 
different kUs for £po/L = 0.5. Similar curves have been 
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generated for other, smaller, values of £po/L. The behav- 
ior of these other curves is the same as the one plotted 
in Fig. |. 

V. CONCLUSIONS 

We conclude with a few brief comments. To fit the 
PE and WLC distributions we have matched the max- 
ima of the two. A fit can also be effected by matching 
moments, or by a least-squares procedure. Alternate fit- 
ting approaches have been explored by us, but the re- 
sults obtained were entirely consistent with the conclu- 
sions set forth above. The reason for the present choice 
was the strong correlation between Odijk's prediction and 
the fit of the distributions. If we match the maxima of 
the two distributions, the electrostatic persistence length 
departs from Odijk's predictions when the distributions 
cease to closely resemble each other, and thus the ratio 
of £e/^Odijk is a good measure of the similarity of the two 
distributions. However, if we used a matching of the first 
or second moment to calculate the effective persistence 
length, the two distributions are clearly distinguishable 
from each other well within the regime in which Odijk's 
formula for the persistence length remains accurate. 

Our results indicate that the difference between the 
radial distribution of the PE and the WLC can be at- 
tributed, at least in part, to the influence of end effects. 
In fact, we believe that the behavior of the persistence 
length is substantially controlled by end effects. One way 
of understanding this is in terms of Odijk's derivation of 
the expression (|^) for the electrostatic persistence length 
1^ . This derivation is based on a calculation of the energy 
of a bent segment of a charged rod. A key assumption 



in this derivation is that the segment takes the form of 
an arc of a circle. End effects are readily associated with 
the difference between the shape of a real bent rod and 
the circular arc assumed in Odijk's derivation. An explo- 
ration of these effects in this context will be described in 
forthcoming work [p^ . 

As noted above, the effect of counterion condensation 
has been ignored throughout the above work. It has been 
shown that counterion condensation modifies the bend- 
ing rigidity of a semiflexible chainjlj, ^ and may result 
in the collapse of the PE chain [|14[. We have performed 
a calculation of the distribution function taking into ac- 
count the attractive interaction due to counterion fluc- 
tuations and observed the signature of collapse. We are 
currently investigating these effects pl|. 

In summary, we have been able to verify that the modi- 
fied WLC as a model of a rod-like PE works in the regime 
in which electrostatic effects play a perturbative role, and 
that the Odijk and OSF forms for the electrostatic per- 
sistence length are quantitatively accurate in this regime. 
We have also found that the above model and formulas 
are accurate in non-perturbative regimes if Debye screen- 
ing is sufficiently strong. In the regimes in which the 
above model fails, we have discovered evidence for a scal- 
ing form, exhibited in Eq. (p^, for the correction to the 
Odijk result for the persistence length. As yet, no model 
provides a theoretical basis for the scaling exhibited by 
this formula. 
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Science Foundation under Grant No. CHE99-88651. 
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FIG. 1: End-to-end distributions of charged semiflexible 
chains. The thin sohd lines are plots of G{r) for PE's with 
kL = 12,6,3 from left to right. In all cases PL = 100 and 
ipo/L = 0.5. The dashed line is the end-to-end distribution 
of a neutral semiflexible chain with Ipo/L = 0.5. 




FIG. 2: The distribution for a polyelectrolyte with Ipo = 0.5, 
/3L = 2400, and kL = 50, compared to that of a wormlike 
chain with Ip = 0.74 (two superimposed curves and perfect 
match with Eq. (0)). Inset: the distribution for a polyelec- 
trolyte with ipd = 0.5, /3L = 5, and kL = 5, compared to 
that of a wormlike chain with fp = 0.524 (two superimposed 
curves and perfect match with Eq. (0)). 



G(r) 




6.6 0.7 0.8 0.9 



FIG. 3: The distribution for a polyelectrolyte with Ipo = 0.5, 
PL — 600, and kL = 10, compared to that of a wormlike chain 
with £p = 1.3. Inset: the distribution for a polyelectrolyte 
with £po — 0.5, /3L = 60, and kL = 10, compared to that of a 
wormlike chain with ip = 0.6 (two superimposed curves). 
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FIG. 4: The diagram delineating the quality of a WLC fit 
for Ipo/L = 0.5. The curve separating the two regimes corre- 
sponds to fe/^odijk = 0.58. This line also indicates the values 
of PL in Eq. at different kL's. 
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FIG. 5: Comparison of our results for the electrostatic persis- 
tence length with Odijk's finite size formula (solid line) 
at kL = 2 with Ipo/L = 0.5 (filled circles) and Ipo/L = 0.01 
(hollow circles). The inset is for nL = 17 and Ipo/L = 0.5. 




FIG. 6: Log-log plot of a comparison of our data, suitably 
rescaled, to the expression in Eq. (^|) The points correspond 
to our rescaled plots of le/L versus f5L at kL = 6, 8, 10, 12. 



